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Ultrashort  Pulse  EflFects  in  Semiconductor  Amplifiers  and  in 

Dispersive  Media 

Arizona  Center  for  Mathematical  Sciences 

R.A.  Indik,  S.VV.  Koch  J.V,  Moloney  and  A.C.  Newell 

The  Arizona  center  for  Mathematical  Sciences  is  currently  involved  in  two  investigations 
that  are  funded  by  AFOSR.  One  of  these  has  been  underway  since  November  1993  and 
involves  the  study  of  ultrashort  laser  pulse  eifects  in  the  eye.  The  other  started  this  February, 
and  involves  study  of  dynamics  in  semiconductor  amplifiers  and  lasers.  Both  of  these  projects 
involve  large  scale  computation,  and  code  development  for  a  massively  parallel  platform. 


Semiconductor  amplifiers  and  lasers 

Our  goal  is  to  use  models  for  the  optical  behavior  of  semiconductors  based  on  an  understand¬ 
ing  of  the  underlying  physics  to  develop  and  test  a  heirarchy  of  eq^uations  and  numerical  codes 
for  modeling  semiconductor  behavior.  The  goal  is  reasonable  because  the  basic  physics  de¬ 
scribing  the  optical  properties  of  semiconductors  is  well  developed.  It  is  challenging  because 
that  physics  happens  over  an  extreme  range  of  time  scales,  from  carrier-carrier  collision  ef¬ 
fects  at  a  femtosecond  scale  to  heat  diffusion  on  a  millisecond  scale.  To  predict  the  behavior 
of  a  device  we  would  like  to  have  model  equations,  derived  from  the  underlying  physics, 
which  do  not  require  that  the  full  range  of  time  scales  be  resolved.  Different  applications 
or  questions  will  require  the  use  of  different  model  equations,  but  all  must  be  mutually 
consistent. 

Our  starting  point  in  this  study  has  been  Maxwell’s  equations  coupled  to  the  semiconduc¬ 
tor  Bloch  equations  (SCB).  The  SCB  equations  describe  the  many  body  interactions  among 
the  electrons  and  holes  in  the  semiconductor  due  to  Coulumb  effects  and  optical  interactions. 
We  are  developing  and  testing  a  Maxwell  solver  to  integrate  full  vector  Maxwell  in  one  two 
or  three  dimensions.  At  this  time  that  equation  is  being  integrated  coupled  to  a  simple  two 
oscillator  material  equation.  Our  next  step  with  this  code  is  to  couple  the  SCB  equations 
into  the  code. 

We  have  also  developed  code  for  integrating  envelope  equations  derived  from  scalar 
Maxwell’s  equations  for  a  scalar  field  coupled  to  the  SCB  equations.  This  code  can  handle 
plane  wave  (one  dimensional)  and  a  single  transverse  dimension,  (two  dimensional).  At  this 
time,  the  codes  are  being  used  to  study  the  behavior  of  ultrashort  pulses  propagating  in  the 
nonlinear  medium.  The  case  of  a  semiconductor  amplifier  has  been  our  focus.  We  have  been 
able  to  integrate  the  plane  wave  equations  for  considerably  longer  propagation  distances  than 
had  been  possible  previously  due  to  the  efficiency  of  our  algorithm.  In  so  doing,  we  have 
predicted  novel  effects  such  as  pulse  compression,  the  appearance  of  an  adiabatic  following 
behavior,  and  spontaneous  pulse  breakup  for  very  high  gain  amplifiers  [?].  In  addition  we 
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have  made  a  study  of  the  interactions  tn  such  an  amplifier  of  simultaneous  ultrashort  strong 
pulses  at  different  frequencies  [?].  This  is  relevant  to  the  ability  of  such  a  device  to  be  used 
with  frequency  multiplexing. 

We  have  also  completed  preliminary  studies  of  the  behavior  of  beams  with  transverse 
structure  in  such  amplifiers.  At  this  time  we  are  focusing  out  efforts  on  improving  the  im¬ 
plementation  of  this  code  on  the  CMo  at  AHPCRC  at  Minnesota  Supercomputing  Center. 
This  has  proved  a  challenging  task.  The  nature  of  the  problem  is  well  adapted  to  a  massively 
parallel  distributed  memory  model,  since  a  domain  decomposition  distributing  calculations 
to  processors  depending  on  the  value  of  the  transverse  spatial  variable  leads  to  calcula¬ 
tions  being  performed  almost  entirely  in  local  memory,  with  extremely  little  cross  processor 
communication  required.  (Only  the  diffractive  term  in  the  field  equation  couples  different 
transverse  spatial  locations  —  the  field  represents  less  than  1%  of  the  data  in  memory,  and 
its  propagation  requires  less  than  0.01%  of  the  CPU  demand  of  the  program.)  However, 
our  initial  efforts  to  implement  this  approach  have  taught  us  that  the  programming  model 
provided  by  CM  FORTRAN  leads  to  a  great  deal  of  cross  processor  communication  that  is 
extremely  diScult  to  control.  We  have  had  some  success  in  improving  the  paradlelization  of 
the  code,  and  have  realized  speedups  using  32  processors  (128  vector  units)  of  a  factor  of  25 
over  single  processor  code,  but  the  size  of  our  problem  is  such,  that  we  feel  it  is  worthwhile 
to  achieve  ma.ximum  parallelization,  which  should  be  very  close  to  a  factor  of  120.  This  is 
particularly  important  since  we  plan  to  study  the  behavior  of  these  equations  when  a  cavity 
is  added,  so  that  the  system  lases.  This  will  require  much  longer  caIc*uIations  than  those 
currently  being  made. 


Ultrashort  pulse  propagation  in  the  eye 

Researchers  at  the  Armstrong  Laboratory,  Brooks  AFB  are  involved  in  setting  safety  stan¬ 
dards  for  ultrashort  laser  pulses.  No  reliable  standard  currently  exists  for  picosecond  and 
sub-picosecond  pulses.  For  femtosecond  pulses  it  appears  that  the  very  high  powers  contained 
in  low  energy  pulses  can  induce  a  nonlinear  self-focusing  of  the  pulse  leading  to  plasma  gen¬ 
eration  [?].  We  are  involved  in  modeling  the  propagation  of  ultrashort  pulses  in  dispersive 
media,  such  as  the  fluid  in  the  eye.  These  equations  are  nonlinear,  and  in  one  approximation, 
they  lead  to  the  two  dimensional  nonlinear  Schrodingcr  equation  with  a  focusing  nonlinear¬ 
ity.  This  equation  has  singularities,  corresponding  to  a  collapse  of  the  optical  beam,  after  a 
finite  propagation  distance.  In  fact,  estimates  of  the  values  of  the  various  parameters  in  this 
simple  model  and  the  powers  used  in  the  laboratory  lead  to  a  distance  to  collapse  somewhat 
smaller  than  the  length  of  the  eye.  Thus  one  possible  mechanism  for  the  observed  threshold 
is  beam  collapse.  As  useful  as  those  simple  models  are  in  getting  a  feel  for  the  behavior  of 
the  system,  there  are  very  important  effects  that  arc  omitted.  Wc  arc  now  investigating  the 

behavior  of  the  system  when  some  of  those  effects  are  included.  In  particular,  the  meium 

_ _ _  / 

^Our  experience  ^dth  the  CMS  is  based  on  a  beta  release  of  the  software,  and  may  not  reflect  the  behavior 
of  later  releases 
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through  which  the  pulse  is  propagating  has  dispersion.  This  dispersion  could  arrest  the  col¬ 
lapse,  since  it  tends  to  spread  the  energy  out  along  the  direction,  of  propagation  at  the  same 
time  as  tiie  collapse  is  concentrating  it  radially.  Whether  the  collapse  is  arrested  or  not,  it  is 
useful  tc  be  able  to  say  what  the  intensities  will  be  within  the  eye.  We  have  derived  simpli¬ 
fied  eq^uations  estimating  these  effects,  and  compared  them  to  direct  numerical  simulations  of 
model  equations  including  dispersion,  diffraction  and  self  phase  modulation  [?,  ?,  ?].  If  the 
collapse  proceeds  sufficiently  far,  nonparaxial  effects  may  become  important.  The  Maxwell 
equation  solver  mentioned  in  the  previous  section  will  be  used  to  study  such  effects.  This  is 
numerically  quite  challenging,  since  the  Maxwell  solver  must  resolve  the  optical  wavelength. 
The  three  dimensional  vector  Maxwell  cede  requires  a  great  deal  of  memcry  ao  well  as  very 
significant  amounts  of  CPU.  This  code  has  been  implemented  in  FORTRAN  90  on  both  the 
Cray  C90  at  WES  and  the  CMS  at  AHPCRC.  It  parallelizes  quite  well,  though  as  in  the 
case  of  the  semiconductor  propagation  code,  it  is  not  a  simple  as  one  would  hope  to  achieve 
the  amount  of  parallelization  potentially  present  in  the  algorithm  while  working  with  CM 
FORTRAN.  Another  potentially,  important  effect  that  we  are  modeling  i.s  the  formation  of 
a  plasma.  The  form  of  the  nonlinearity  in  a  plasma  is  defocusing,  so  if  a  plasma  forms  (due 
to  the  high  intensities  of  the  self  focusing  beam),  it  may  arrest  that  collapse. 
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